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Abstract 

A new calculus on fractal curves, such as the von Koch curve, is for- 
mulated. We define a Riemann-like integral along a fractal curve F, called 
_F"-integral, where a is the dimension of F. A derivative along the fractal 
curve called _F"-derivative, is also defined. The mass function, a measure- 
like algorithmic quantity on the curves, plays a central role in the formula- 
tion. An appropriate algorithm to calculate the mass function is presented 
to emphasize its algorithmic aspect. 

Several aspects of this calculus retain much of the simplicity of ordi- 
nary calculus. We establish a conjugacy between this calculus and ordi- 
nary calculus on the real line. The _F"-integral and F"- derivative are 
shown to be conjugate to the Riemann integral and ordinary derivative 
respectively. In fact, they can thus be evalutated using the corresponding 
operators in ordinary calculus and conjugacy. Sobolev Spaces are con- 
structed on F, and _F"- differentiability is generalized . Finally we touch 
upon an example of absorption along fractal paths, to illustrate the utility 
of the framework in model making. 



1 Introduction 

It is now well known that fractals pervade nature [U [5] . The geometry of fractals 
is also well studied [2 |31 HI [SJ 111 [7] . Fractal curves often lack the smoothness 
properties required by ordinary calculus. For example, observed path of a quan- 
tum mechanical particle [8j or Brownian and Fractional Brownian trajectories 
[U [3] are known to be fractals and are continuous but non-differentiable. A 
percolating path, just above the percolating phase transition can be considered 
as an appoximate realization of a fractal curve [3] . If a long polymer is modeled 
as a fractal curve, then accumulation of a physical property along the curve 
would amount to integration on such a curve. This is often carried out using ad 
hoc procedures. 

While there are some remarkable approaches to develop tools for such situa- 
tions [ini El HH [131 m] , much more is desired. This paper aims to formulate a 
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calculus specifically tailored for fractal curves, in a close analogy with ordinary 
calculus. In particular, we adopt a Riemann-Stieltjes like approach for defining 
integrals, because of its simplicity and advantage from algorithmic point of view. 
Such an approach was concieved in (15) which began with formulation in the 
terms of Local Fractional derivatives. A new prescription was proposed to give 
meaning to differential equations on Cantor-like sets which are totally discon- 
nected where the Local Fractional Derivatives do not carry over. The calculus 
formulated and developed in [1^1 [HI [TH] for fractal subsets of R fully justifies 
the prescription in |15| . In particular, an integral and a derivative of order a are 
defined [TB] on Cantor-like totally disconnected subsets of the real line, where 
a 6 (0, 1] is the dimension of F. This calculus, called F"- calculus has many 
results analogous to ordinary calculus and can be viewed as a generalization of 
ordinary calculus on R. In fact, in [T7l[l8] a conjugacy between the i^"-calculus 
and ordinary calculus is discussed. 

The present paper extends that approach, which was developed for discon- 
nected sets like Cantor-sets, to formulate calculus on fractal curves which are 
continuous . The organization of the paper is as follows. In Section[5]we define a 
mass function and integral staircase function. The mass function gives the con- 
tent of a continuous piece of the fractal curve F. The staircase function, more 
appropriately called the "rise function", is obtained from the mass function and 
describes the rise of the mass of the curve with respect to the parameter. We 
emphasize the algorithmic nature of the mass function: by presenting an algo- 
rithm to calculate it. In section [3] we show that the mass function allows us 
to define a new dimension called j-dimension, which is algorithmic and finer 
than the box dimension. In section |4] we discuss the algorithmic nature of mass 
function and present an algorithm to calculate it. In section [STTl the concepts of 
limits and continuity are adapted to the concepts of F-limit and F-continuity. 
Section 15.21 is devoted to the discussion of integral on fractal curves called F"- 
integral. The formulation is analogous to the Riemann integration |19| . The 
notion of F"-differentiation is introduced in section 15.31 The fundamental the- 
orems of F"-calculus proved in section [531 state that the ^"-integral and F"- 
derivative are inverses of each other. The conjugacy between ^"-calculus on 
F and ordinary calculus on the real line, discussed in section [HI establishes a 
relation between the two and gives a simple method to evaluate _F"-integrals 
and F"- derivatives of functions on the fractal F. In section function spaces 
of F"-integrable and F"-differentiable functions on the fractal F are explored. 
In particular Sobolev Spaces are introduced and abstract Sobolev derivatives 
are constructed. Finally as a simple physical application we briefly touch upon, 
as an example, a simple model of absorption along a fractal path in section [71 
Section [5] is the concluding section. 

2 The mass function and the staircase 

This paper can be considered as a logical extension of calculus on fractal subsets 
of the real line developed in [18]. The proofs which are analogous to those in 
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[15] are omitted. 

In this paper we consider fractal curves, i.e. images of continuous functions 
/ : R — > R" which are fractals. To be precise: 
Let [qq, bo] be a closed interval of the real line. 

Definition 1 A fractal (curve) F C R" is said to be continuously parametriz- 
able (or just parametrizable for brevity) if there exists a function w : [ao, bo] — >■ 
F C R" which is continuous, one-to-one and onto F. 

In this paper F will always denote such a fractal curve. 
Examples: 

1. A simple example of such a parametrization is the function w : R ^ R^ 
defined by w(t) = {t,W^{t)) where W^{t) is the well known Weierstrass 
function p] given by 

oo 

W^{t) ^^X^"-^')'' sin XH 

k=l 

where A > 1 and 1 < s < 2. The graph of W^{t) is known to be a fractal 
curve with box-dimension s. 

2. Our next example constitutes of one important class of parametrizations of 
self-similar curves in two dimensions (There are other ways of parametriz- 
ing fractal curves ; for example see [201)- Let T^, « = 0, . . . , n — 1 be linear 
operations which are composed of rotation and scaling. Each Ti can be 
represented by a 2 x 2 matrix: 

rp _ g cos 9i — sin 9i 
' * sin^i cos^i 

Further, they should satisfy the condition: 

n-l 
1=0 

for any vector v, and < s,; < 1 for i = 0, . . . , n — 1. The fractal is defined 
by the limit set [3] of the similarity transformations: 

i-i 

SjM = + J - 0, . . . , n - 1 

i=0 

where Vq is a fixed vector. The limit set will be in the form of a curve 
because of the way Sj are constructed from T^. 

Let lnt\ denote the integer part of nt. Now, the function w defined 
implicitly by 

w(t) = J2 r,(vo) -|-TL„tj(w(nt - lnt\)), 0<t<l (1) 

2 = 
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parametrizes the above fractal curve. To implement it as an algorithm, 
we stop the recursion at some appropriate depth. The continuity and 
invertibility of this parametrization can be numerically verified, when the 
curve itself is non-self-intersecting. 

In particular the von Koch curve is realized by setting all Sj = 1/3, = 
6*3 = 0, 61 = —62 = 7r/3, and vq = (1, 0) (the unit vector along x axis). 



Hereafter symbols such as a,6,c,etc denote numbers in [qq, bo] and 9, 0' etc denote 
points of F. 

Definition 2 For a set F and a subdivision P[a,fe],a < b, [a,b] C [aoj&o] 

^ ' ^ ^ T(a + 1) ^ ^ 

i=0 \ < J 

where | • | denotes the euclidean norm on R", 1 < a < n and P[a,b] = {a = 
to, ■ ■ ■ ,tn = b}. 

Next we define the coarsed grained mass function. 

Definition 3 Given 6 > and ao < a < b < bo, the coarse grained mass 
^^{F,a,b) is given by 

^t{F,a,b)= iiif a"[F,P] (3) 

|-P[a,i>]:|-P|<0} 

where \P\= max {ti+\ — ti) for a subdivision P. 

0<i<n—l 

The mass function is the limit of the coarse-grained mass as (5 : 
Definition 4 For ao < a < b < bo, the mass function ^"{F, a, b) is given by 

r{F,a,b) = \imjnF,a,b) 

(5->0 

Remark: Since 7 is a monotonic function of 6. The limit exists , but could 
be finite or +00. 

The following properties of the mass function follow easily. 

Properties of 7"(F, a, b) 

• For ao<a<b<c<bo and 7" {F, a,c) < 00 

7"(F,a,c)=7"(F,a,6)+7"(i^,&,c). 

• 7"(i^, a, b) is increasing in b and decreasing in a. 
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• If 7"(F, a, b) is finite, j°'{F, a, t) is continuous for t £ [a, b]. 

Remark: The implication of this result is that no single point has a nonzero 
mass, or in other words, the mass function is atomless. 

• Let F C R" be parametrizable. Let A be a positive real number, v e R", 
and let T be a rotation operator. We denote 

F + v = {w{t) + v:t€ [ao,bo]} 
XF = {Xw{t) :tG [ao,bo]}. 

and 

TF={Tw{t):t€[ao,bo]} 

Then, 

1. Translation : 



2. Scaling : 

3. Rotation : 



riF + ^r,a,b)=r{F,a,b) 
r{^F,a,b) = X"nF,a,b) 
r{TF,a,b)=r{F,a,b) 



Re-parametrization Invariance of Mass Function 

The definitions of a", 7^, and therefore 7" implicitly involve the particular 
parametrization w. Here we show that although defined through the parametriza- 
tion, these definitions arc invariant under the change of parametrization. In or- 
der to be able to unambiguously and explicitly refer to the parametrization, we 
introduce a temporary change in the notation to explicitly indicate dependence 
on parametrization. Thus given a parametrization w : [a,b] — >■ R", we use the 
following notation here: 



a"[F,P;w] 



!^^ |w(f,-+,)-w(f,)r 



75"(i^,a,6;w) = inf a«[F,P;w] 

|P|<5 

7"(F,a,6;w) = limjf{F,a,b;v,) 

S->-0 

Let wi and W2 be two parametrizations of the given fractal curve. By our 
definition of parametrization, Wi and W2 are continuous and one-to-one. Let the 
domain of Wi be [ai, 61], and that of W2 be [a2, 62]. We further assume that Wi 
and W2 have the same orientation, i. e. wi(ai) = W2(a2) and wi(6i) = W2(62). 
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Thus, z = w^^ o wi : [fli, 61] — > [a2, 62] is a continuous, one-to-one and strictly 
monotonically increasing function. 

Now, given i52 > and e > 0, there exists a subdivision P2 of [02,62] such 
that 

cr" {F, P2 ; W2] < "2 , ^2 ; W2 ) + e. 

The set of points P\ = {z~^{t) : t G P2} forms a subdivision of [ai, 61]. Then, 

a«[F,Pi;wi] =a"[F,P2;w2] 

by appropriate substitution. Therefore, 

a"[P, Pi; wi] < (P, 02, 62; W2) + e 

which impUes that 

7"^(P,ai,&i;wi) < 7fc(-^'a2,62;w2) +e 

where = |Pi|. Further, since z is continuous, hm(5i = 0, as 82 — ^ 0, implying 
that 

7"(P,ai,6i;wi) < 7"(P, 02, 62; W2) + e. 

Since e is arbitrary, and the same argument remains valid starting with z~^ = 
wj"^ o W2, we conclude that 

7"(P,ai,6i;wi) = 7'*(P,a2,62;w2). 

This establishes the fact that the mass function depends only on the fractal curve 
(i. e. the image of the parametrization), and is independent of the parametriza- 
tion itself. Since the mass function underlies the calculus developed in the subse- 
quent sections, the calculus is also independent of the particular parametrization 
chosen. 

Now we introduce the integral staircase function for a set F of order a. 

Definition 5 Let po S [aoi^o] be arbitrary but fixed. The staircase function 
Sp : [ao, 60] -^H of order a for a set F is given by 

cccf.-. _ / T{F,Po,t) t>po f.. 
j -jC.^F,t,po) t<po 

where t G [ao, 6o]- 

In the rest of this paper we take po = ao unless stated otherwise. 

Here this function may, more appropriately, be described as a rise function. 
However we retain the name staircase function because in analogous calculus on 
fractal subsets of the real line this role is played by a staircase. 

Throughout the paper we consider only those sets for which Sp is strictly 
increasing and thus invertible. Further, we define 

J{0)=SUw-'{0)), e&F (5) 
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staircase 




Figure 1: Sp for von Koch curve. The von Koch curve lies in the XY plane. 
The vertical lines are drawn to guide the eye (to show how Sp rises) 

which is the function induced by Sp on F, and it is also one-to-one. 

As an example, figure [1] shows the staircase function for the von koch curve. 
The curve was parametrized as given in |20| . 

A log- log graph of the staircase function Sp{t) against the Euclidean distance 
between origin and w(t) for the von- koch curve is shown in fig [2] 

3 The 7- Dimension 

We now consider the sets F for which the mass function ^"{F,a,b) gives the 
most useful information. Due to the similarity of the definitions of mass function 
and the HausdorfF outer measure, the former can be used to define a fractal 
dimension as follows. 

It can be seen that ^"{F, a, b) is infinite upto certain value of a, say ag, and 
jumps down to zero for a > ao. Thus 

Definition 6 The 'y- dimension of F, denoted by dim-y(i^), is 

dim..,(F) = inf{a : -f°'{F,a,b) = 0} = sup{a : j"{F,a,b) = 00} 
It follows that the 7-dimension is finer than the box dimension. Thus 

dim^{F) < dims(i^). 
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Figure 2: log-log graph of Euclidean distance between origin and w(t) (Y-axiz) 
vs Sp{t) for i G [0, 1] for von-koch curve 

7-dimension for self-similar curves 

Let a denote the 7- dimension of a self similar curve , which is made up of m 
copies of itself, scaled by a factor of ^ and rotated and translated appropriately. 
Then using the translation, scaling and rotation properties of the mass function, 
one can see that the mass of the whole curve is given by 



7"(F,ao,6o) = m7"(-F,ao,&o) 
n 

r{F,ao,bo) ^ m{-mF,ao,bo) 



(6) 




(7) 



This is same as the Hausdorff dimension of self-similar curves [4] . 
Thus for self-similar curves 



dim^ F = dim-H F = dim b F 



where dim-^F denotes the Hausdorff dimension and dinriBF the box dimension 
of F. 



4 Algorithmic Nature of the Mass Function 

Let us first summarize the definition of mass function: 



One of the main difference between the Hausdorff measure and the mass function 
is that while the HausdorfF measiire is based on sums over a countable covers 
(composed of arbitrary sets) of the given set F, the mass function is based on 
finite subdivisions of the parametrization domain. Prom an algorithmic point 
of view, the extent of the set of all possible finite subdivisions is much smaller 
than that of all countable (finite and infinite) covers of a set. This makes the 
mass function much more amenable to an algorithmic computation. 

As in any algorithm which intends to approximate the infimum, we would 
like to find a subdivision P such that cr°'[F, P] is close to the infimum. Further, 
we can consider values of 5 only as small as practically possible within the reach 
of numerical calculations. The goal of the algorithm is thus to find a subdivision 
P as described above, given a fixed 5. 

However, the set of allowed subdivisions is still large, to explore all of it 
systematically. Further the constraint |P| < 6 does not restrict the number 
of points in P, rendering the standard deterministic optimization algorithms 
either inapplicable or too complex to implement. More appropriate is a Monte 
Carlo method where a subdivision is modified in a variety of ways randomly but 
consistently with the constraint |P| < 5, and the change is accepted if the sum 
(t" [F, P] decreases due to the modification. The algorithm presented below, is 
based on this strategy. 

A Monte Carlo Algorithm 

For the purpose of this algorithm, [a, b] denotes tlie domain of w. Further, 
"randomly" means with a uniform probability unless stated otherwise. The 
symbol P always indicates the "current" subdivision in consideration. 

We begin with a uniform subdivision P such that |P| = (5/4, and iteratively 
improve it using the following prescription. 

1. Choose two numbers x,y € [a,b] randomly, and relabel them if necessary 
so that X < y. Then [x,y] C [a,b]. Let P' ~ {ti : < i < m} denote the 
set of all points of P fl [x,y]. We now modify P' in one of the following 
ways with equal probability, and denote the resultant by P": 

(a) With a probability Pc = min(l,(5/(t/ — a:)), we shift each point U 

(except to and tm) by a random amount between [—6/2,6/2], if the 
resultant subdivision P" still satisfies \P"\ < 6. 

(b) With a probability pd = min(l,(5/(y — x)), we remove each point 
ti (except to and tm) from P', if the resultant subdivision P" still 
satisfies |P"| < 6. 
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Figure 3: The image (under w) of a numerically computed near-optimal subdi- 
vision P, for 6 = 0.05, superimposed on the von-Koch curve. 

(c) With a probability pi = min(l, 5/{y — x)), we insert a point between 
each ti and U+i which is chosen randomly from [ti,ti+i]. (However, 
to avoid accumulating too much of rounding error, we insert the point 
only if the distance between ti and U+i is greater than 6/10.) 

2. Form a new subdivision Pi = {P n [a,x)) U P" U (P n {y,b]), i. e. the 
subdivision of which the points belonging to [x,y] are changed by the 
above procedure. If a"[F,Pi] < cr"[P, P], then we consider Pi as the 
"current" subdivision which will be possibly improved further using above 
steps. Otherwise we consider P again for the purpose. 

As the sum cr" [F, P] approaches the infimum, many of the newly formed 
subdivisions P' are rejected since they sum up higher than P. Thus near the 
infimum, the sum remains constant for many consecutive iterations, and changes 
only intermittently. Therefore the usual convergence criterion of terminating 
iteration when the difference between successive iterations or every K iterations 
{K being a suitable large integer) goes below certain small number, is not useful 
in this case. Instead, after examining the sum over a large number of iterations, 
we observe that the sum stops making significant progress between N' — 1000 
to N' = 2000, where TV' = N/n is the number of iterations N normalized by the 
current subdivision size n. Further, we need to go through all these iterations 
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Figure 4: The evolution of r{a + l).a°'[F, P] over the normalized number of 
iterations. The evolution is shown only up to N' — 100, since the latter part 
(100 < A^' < 2000) is almost flat and uninteresting. 



more than once, just to ensure that subdivision is really optimal. Occasionally 
it may happen that the sum settles a little above the optimal value, gettting 
"trapped" in a "local minimum". 

We demonstrate the results of this algorithm as applied on the von Koch 
curve, parametrized as in equation ([T|). It turns out that the mass of the entire 
von Koch curve is a little less than 0.51/r(Q; + 1), a = ln(4)/ln(3). The image 
(under w) of the optimal subdivision found by the algorithm is shown in fig- 
ure [31 superimposed on the von Koch curve. The evolution of the sum over the 
normalized number of iterations is shown in figure HI 

The above description assumes that the value of a is the same as the 7- 
dimension of the set F, say aQ. We expect (5- independence in the values of 
(7"[F, P((5)] where P{S) denotes the resultant subdivision of the algorithm at 
the scale S, since the value of 7^ converges to a finite nonzero value. This is 
what we observe from the values of (t"[F, P{5)] obtained for various values of 5 
(figure ID). 

Now we would like to consider cases when a ^ ag. Let < Si < 62- If 
a < ao, then j°'{F,a,b) — 00. Therefore we expect that R{a) = a^lF, P{6i)] 
/a^lF, P{S2)] > 1. Similarly since a > ao implies "f°'{F,a,h) — 0, we expect 
that R{a) < 1. 
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This fact can be used to algorithmically calculate the 7-dimension ag: We 
need to find the number ao such that R{ao) = 1. We already know that 
ao e [l,m], m being the embedding dimension, since F e is a curve. 
Treating this as the initial bracket of values for ao, wc just need to use some 
algorithm such as bisection to shrink this bracket to sufficient accuracy. 

5 The F"-Calculus 

Most of the proofs which are similar to the proofs in the case of discontinuous 
sets like Cantor-like sets are omitted. 

5.1 F-Limit and F-Continuity 

Now we introduce limits and continuity along a fractal curve. 

Definition 7 Let F C R" be a fractal curve, and let f : F ^ R. Let 9 e F. A 
number I is said to be the limit of f throught points of F, or simply F- limit, 
as 6' — >■ 9, if given e > there exists 5 > such that 

9' eF and \e' - e\ < 6 ^ 1/(6*') - l\ < e 

If such a number exists it is denoted by 

Z = F-lim^/(^') 

Definition 8 A function f : F ^ R is said to fee F- continuous at 6 £ F if 
f{9) = F -lime- ^ef{0'). 

Definition 9 f : F ^ R is said to be uniformly continuous on E C F if for 
any e > there exists 6 > such that for any 9 £ F and 9' £ E 

\e'-9\<S^\f{9')-f{9)\<e 

5.2 F°-Integration 

Wc denote the class of bounded functions f : F ^ Rhj B{F). 

Definition 10 For ti,t2 € [ao,ho\,ti < t2 a section or segment C{ti,t2) of the 
curve is defined as 

C{ti,t2) = {vf{t'):t' G[ti,t2]} 
Definition 11 Let / : F — >■ R and t\,t2 € [ao,6o];ii < ^2 o-nd let 
M[f,C{tut2)]= sup f{9) 

6eC(tiM) 

and 

m[f,C{tut2)\= ^ inf f{9) 
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Let Sp(t) be finite for t G [a, b] C [ao, 60] • Let P — {to, ■ • • , tn} be a subdi- 
vision of [a,b] . The upper and the lower F^-sum for the function f over the 
subdivision P are given respectively by 



n-l 



U"[f,F,P] = Y,M[f,C{t,,t,+,)][SUU+i) - S^{t,)], (9) 



i=0 
n-l 



L"[/, F,P] = J2 Hf: Cit,, t,+i)][S^{t,+i) - S^{t,)]. (10) 

1=0 

Now we define the F"-integral 

Definition 12 Let F be such that Sp is finite on [a, b]. For f G B{F), the lower 
and upper F" -integral of the function f respectively , on the section C{a, b) are 

[ f{e)d''pe^ sup L"[f,F,p] (11) 

JC(a,b) P[c,,b] 

[ f{e)d"pd= vai U''[f,F,P] (12) 
If f G B{f), we say that f is F'^- integrable on C{a,b) if 



f{6)d^p6= / f{e)d%d 

C{a,b) JC{afi) 

and the common value is called the F" -integral, denoted by 

f{e)d%e. 



C{a,b) 



5.3 F°-Differentiation 



Definition 13 Let F be a fractal curve. Then the F" -derivative of function f 
at 6 E F is defined as 



if the limit exists. 



Theorem 14 If{Dpf){6) exists for all 9 G C{a,b), then f is F-continuous on 
C{a,b). 

Remark: The F"-derivative Dp{f) of a constant function f : F ^ R, 
f{9) = k G Ris zero. This result is to be contrasted with the classical fractional 
derivative (Riemann-Liouville, and others) of a constant, which is not zero in 
general [H HI [H [H] . 
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5.4 Fundamental theorems of F'^-calculus 

The F"-integratioii and F"-difFerentiation are related as inverse operations of 
each other. The first fundamental theorem states: 

Theorem 15 Let f G B{F) is an F-continuous function on C{a,h),and let 
g : f ^ H be defined as 



The second fundamental theorem says that the i^"-integral as a function of 
upper limit is the inverse of .F^-derivative except for an additive constant. 

Theorem 16 Let f : F ^ R be continuously F"- differ entiable function and 
h : F be F-continuous, such that h{6) — {Dpf){0)). Then 



Comparison with approaches involving Local Fractional Op- 
erators 

The local fractional derivative (LFD) operator constructed in 25 was based on 
the renormalization of Riemann-Liouville differential operator on the real line. 
It was utilized to establish the relation between the differentiability properties 
of nowhere differentiablc functions, such as Weierstrass function, and the di- 
mension (Holder exponent) of its graph. The domain of these functions is R 
and not a fractal. In [THj though the prescription was developed using LFD, 
it was realized that to make the local fractional Fokker Planck equation causal 
and dynamically consistent, the evolution had to be restricted to fractal subsets 
of the real line. Moreover the order of differentiation had to be the dimension 
of the fractal support. Since Cantor- like sets are totally disconnected, this 
necessitated a rigorous development of calculus on fractals from first principles, 
without using standard fractional calculus on R, which was carried out in [16] 
The present paper is a logical extension of such a formulation, to fractal curves. 
While Cantor-like sets considered in [12] are totally disconnected, the von-Koch 
like fractal curves are continuous but tangentless. Thus all the constructions 
and proofs have to be carried out keeping in mind this difference of the domain 
of functions and operators. (Thus for example the notion of 'set of change' and 
'a-perfect sets' was crutial for cantor-like sets in [TB] which is repleced by the 
invertibility of Sp in case of fractal curves considered in this paper.) 




for all t e [a, b] . Then 



{D'^Fdm = f{e) 



and 9 — w(t) 
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There is multiplicity of approaches leading to the notion of local fractional 
calculus. Various authors [26], [27], [28] have further developed the notion of 
local fractional differentiation with different approaches. They provide suitable 
framework for different classes of problems. The development in [33] is based 
on approach involving difference operators. In finding derivative the quotient is 
taken with respect to /i" where h is an increment of the independent variable. 
This can be contrasted with the use of Sp in the present approach. This also 
reflects in the Taylor series where powers of (Sp) appear (see equation (fTTjl 
below) rather than powers of h itself as in [53] . Further the domain is R in [53] 
whereas it is a fractal curve in the present paper. 

In [57] the notion of classical fractional derivative is modified. Again the 
essential difference mentioned above for Taylor series and the domain functions 
is also to be noted here. The development in [28] is based on the Weyl Derivative 
and the domain of functions is R and not a fractal curve. 

We may emphasize that in the present approach the role of the independent 
variables is delegated to the staircase/rise function Sp, see e.g equation ([T^ . 
In this sense our approach is like Stieltjes approach in spirit but with some 
essential differences as noted in [16]. The function Sp captures the essence of 
fractal support, hence its use makes the calculus suitable for fractals. 

6 Conjugacy of F"-Calculus and Ordinary Cal- 
culus 

In this section, we define a map (j) which takes an F"-integrable function / : 
F — > R to a Riemann integrable function g : [Sp{ao), Sp{bo)] R such that 
their corresponding integrals have equal values. Thus, the map (j) exhibits a 
conjugacy between the two operations. 

First let us define certain classes of functions: 

1. B{F) : class of bounded functions / : — s> R. 

2. B{[c,d]) : class of bounded functions / : [c, d] — > R 

3. C{F): set of all functions which are i^"-integrable on C(ao, 6o). 

4. The image of F under Sp is denoted by i.e K — [S'^(ao), 'S'J? (6o)], and 
B{K) denotes the class of functions bounded on K . 

5. C{K) denotes the class of functions in B{K) which are Riemann integrable 
over the interval K = [S'^(ao), 5^(60)]- 

In order to fix the notation, here we briefly review the definition of Riemann 
integral. Firstly, if g G B{\c,d\) and / C [c, d] is a closed interval, then we 
denote M'[g,I] = s\vg^^j g{x) and m'[g,I] = vaixei g{x)- Further, the upper 
and lower sum over a subdivision P[c,d] = {yoj ■ • ■ tVu] are given by U'\g,P] = 
YJ^Zo M'igAyi^Vi+i]] and L'{g,P] = ^"=0^ [^ii 2/«+i]]- If the upper and 
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lower integrals given respectively by inip U'[g, P] and swppL'[g,P] are equal, 
then g is said to be Riemann integrable, and the Riemann integral 




is defined to be the common value. 

Now we define the above mentioned map (p: 



Definition 17 The map cf) : B{F) B{[S^{ao), S^{bo)]) takes f G B{F) to 
(j)[f] G B{[Sf{ao), S^{bo)]) such that for each t G [ao, bo], 



Lemma 18 The map <p : B{F) — > B{K) is one to one and onto. 

The proof is straightforward. Thus we are assured that the inverse map (f)~^ 
exists. 

The following theorem brings out the conjugacy between F"- integrals of 
functions along the fractal curve F and the Riemann integrals of their images 

under (j). 

Theorem 19 A function f G B{F) is F" -integrable over C{a,b) if and only if 
g = </)[/] is Riemann integrable over [Sp{a),Sp{b)]. In other words, a function 
f G B{F) belongs to jC{F) if and only if g € ^{K). Further 



Proof: Let f : F ^ R he F"-integrable. Then there exists a subdivision 
P[a,b] = {to,ti,...,tn} such that 



for any e > 0. 

Denote pi = Sp{ti). Then Q = {y, : < i < n} is a subdivision of 




U^[f,F,P]-L^[f,F,P]<e 



(14) 




M[f,C{ti,ti+i)] = 



sup /(w) 



wec(ti,ti+i) 



sup /(w(i)) 



= sup g{S^{t)) 



t6[t»,t, + i] 



sup g{y) 



y^lyi,yi+i] 



= M'[g,[yi,yi+i]] 
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Therefore, 

n-1 

U^[f,F,P] = ^M[/,C(t„t,+i)][5^(i,+i)-5^(t,)] 

1=0 
n-1 

= ^M[f,C{U,U+i)][y,+i -y^] 



i=0 

n-1 



i=0 

= U'[g,Q] (15) 

Similarly 

L"[/,F,P] =L'[5,g] (16) 
then using equations ([14]), (fTS)) and (IT6l) 

f/'[5,g]-i'[.9,g]<e 
which implies that (/ is Riemann integrable over [Sp{a), Sp{b)] and 

g{u)du = / /(0)d^0 

Conversely if 5 is Riemann Integrable, then for given e > there exists a sub- 
division Q' = {vo,...,Vm} of [S'^{a),Sf{b)] such that U'[g,Q'] - L'[g,Q'] < e. 
Then the converse can be proved by following the above steps in the reverse 
order. 
• 

Let /i denote the indefinite F"-integral viz. /i(w(t)) = J^^^ f{9)dp0 and 
let gi denote the ordinary indefinite Riemann integral viz. gi {y) = J^^ g{y')dy'. 
If we further denote the indefinite _F"-integral operator by Ip and the indef- 
inite Riemann integral operator by /, then the result of theorem (fT9)) can be 
expressed as 

as displayed in the commutative diagram of figure El 

The following theorem brings out the conjugacy between i^"-derivative and 
ordinary derivative. 

Theorem 20 Let h be a function in B{F) such that g = is ordinarily 
differentiable on K = range of Sp. Then Dph{9) exists for all 9 E F and 

Upll(p) — \v = J{9) 

Proof: Let v G K. Then by definition 

dg ,. g{u) - g{v) 

— = lim 

dv u-^v u — V 
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i.e given eo > 0, there exists do > such that 

\u-v\<6o^\- < eo 

dv u - V 

Let us recall our assumption that Sp is monotonically increasing and one-to- 
one. Let t = {S<^)-^{v), t' = (S'^)-i(u). Then t,t' G [ao, 6o],/i(w(f )) = g{u) 
and /i(w(t)) = g{v). Thus, 

Since (w)~^ and Sp are continuous, so is their composition Sp o (w)~^. There- 
fore, there exists (5i > such that 

|w(i') - w(i)| <6i^ \S'^{t') - S^{t)\ < 6o 

dg h{w{t') - h{-w{t)) 



'dv S^{t')-S^it) ' 
Setting 9' = w{t'),9 = w{t), we can rewrite this as 



< eo. 



which by definition of f -limit and Dp means 



Theorem 21 Let h G B{F) be an F" -differentiable function at all 9 £ F. 
Further, let g = (l)[h]. Then dg/dv exists at v = J{6) and 

-^\v=j{e) - ■L>ph(9) 

Proof: Asg = cplh], we have g{Sf{t)) = /i(w(t)) for all t G [ao,6o] i-e g{J{0)) = 
h{9) for all 9 eF. 

By definition and substitution 

D^phi9) = F-lim 

^ ^ ' J(9') - J{e) 

J{9') - J{9) 

Thus given eo > there exists ^o' > such that 

\9' 91-,/^'^ g{J{e'))-g[J[0)) . ^ 
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Let V = J {9) and u = J [9'), i.e 9 = J-^v) and 9' = J-^{u). Since J'^ is 
continuous, there exists a (5 > such that 



\u-v\<5=^\9' -9\< 5o 



g{u)-g{v) 



u — V 



D'^h[9)\ < Co 



Which by definition of ordinary derivative gives 



dg 



= lini 9{n)-9{v) ^ ^^^^^^^ 



u^v U — V 



This conjugacy can also be expressed as = 
commutative diagram of figure [5] 



'^D4> as shown in the 



91 



51 



D 



Figure 5: The relation between i^"-integral and Riemann integral, also between 
F"-derivative and Ordinary derivative 



6.1 Taylor Series 

One can write a fractal Taylor series for functions on fractal curve F, by using 
the results of this section. 

If g = 4>[h] be such that the ordinary Taylor series is given by 

ai-) = L 

n=0 

is valid for w, y G [Sp{a),Sp{h)], then for 9,9' & F it can be seen that 

provided h G B{F) is F"- differentiable any number of times on C{a,b) such 
that {D'^y'h e B{F) for any integer n > 0. 
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6.2 Function Spaces in F°-Calculus 

We introduce the foUowing spaces: 

The set of aU functions that have F-continuous i^"-derivatives upto order k 
can be analysed analogous to [T7] and are defined by C''{F),k e N:Set of all 
functions / : F — ^ R such that 

{D'^Tf e C°(F) for al\n<k 

One can define norm on C^{F) for F C R", similar to what is defined for F C R 
using the F"-derivative as follows: 

11/11= ^^v\mTm f&c\F) 

C^{F) are complete with repect to this norm. Unlike in [IT the need for Sp- 
concordant functions does not arise since Sp{u) ^ Sp{v) even if f{w{u)) = 
f{w{v)) for u ^ V. It can also be shown that quite easily that C^{F) is separa- 
ble. 

The spaces of F"-integrable functions and their completion can also be con- 
structed in an analogous manner as is done in ^I6i for fractal subsets of real 
line. 

Set C{F) of F"-integrable functions is a vector space with usual operations 
of addition and scalar multiplication. An appropriate norm Mp can be defined 
for F" integrable functions which satisfies all the required properties. 

AAp(/) = ||/||p = [/ \f{e)\pd%eY'p i<p<oo 

JC(a,b) 

Ap can be shown to act as a norm on L'{F) where L'{F) is a vector space 
of equivalent classes of C{F), the class of all F"- integrable functions. L'^{F) 
which is F'(F) with specific norm Mp is not complete but can be completed 
using standard procedure. The complete space is then denoted by Lp{F) and is 
a Banach space. The spaces Lp{F) and Lp{F) can also be shown to be separable. 

Analogues of abstract Sobolev Spaces can be constructed in exactly the same 
way as is done in the above cited reference for subsets of Real line. 

7 Example: Absorption on fractal curves 

Consider the flux of a fluid or of particles moving steadily through and getting 
absorbed in a percolating cluster or fractured rock. A simple mathematical 
model of this process for a single branch would be that of particles getting 
absorbed along a fractal path. The absorption process can then be modelled by 
the following equation: 

D'^pPie) = ~^p{e) (18) 

where p{0) is the density of fluid at a point of the fractal channel (e.g. backbone 
of the percolating cluster), n being the coefficient of absorption (which in a 
simple model is taken as constant), and D'^ is the F" - derivative. 
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The left hand side of equation ([THl) represents the space rate of change of 
density of particles at a particular position on the fractal curve (or path) varying 
along the path. 

The exact solution of the above equation can be obtained using the conjugacy 
between ^""-derivative and ordinary derivative and applying the corresponding 
operator on p{9) as follows: 

(j)p{0) — p{y) where y — J{9) 

then equation becomes 

ay 

the solution of which is given by 

P{y) = P(0)exp(-Ky) 

Applying the inverse conjugate operator to the above equation we obtain 

p(0)=p(O)exp(-AtJ(^)) 

or since 6 ~ w{u) and J{0) — Sp{u), 

p(w(u)) =p(0)cxp(-At5^(u)) (19) 

This is like stretched exponential behaviour in view of the relation between 
euclidean distance and staircase (see fig [2]). 

8 Conclusion 

In this paper we have developed a calculus on parametrizable fractal curves of di- 
mension a £ [1, n]. This involves the identification of the important role played 
by the mass function and the corresponding (rise) staircase function which may 
be compared with the role played by the independent variable itself in ordinary 
calculus. The definitions of F"-integral and i^"-derivative are specifically tai- 
lored for fractal curves of dimension a. Further they reduce to Riemann integral 
and ordinary derivative respectively, when F — R and a — 1. 

Much of the development of this calculus is carried in analogy with the 
ordinary calculus. Specifically, we have adopted Riemann-Stieltjes like approach 
for integration, as it is direct, simple and advantageous from algorithmic point 
of view. The example of absorption on fractal curves mentioned in section [7] 
demonstrates the utility of such a framework in modelling. Other applications 
may include fractal Langevin equation for Brownian motion and Levy processes 
on such curves, which will follow in future work. This approach may be further 
useful in dealing with path integrals and other similar applications. Another 
direction for extension of the considerations in this paper is the extension to 
crumpled or fractal surfaces which are continuously parametrizable by a finite 
number of variables. 
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